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Abstract 

If a many-body quantum system approaches thermal equilibrium from 
a generic initial state, then the expectation value {il){t)\Ai\ip{t)) ^ where 
\'4){t)) is the system's state vector and Ai is an experimentally accessible 
observable, should approach a constant value which is independent of 
the initial state, and equal to a thermal average of at an appropriate 
temperature. We show that this is the case for all simple observables 
whenever the system is classically chaotic. 
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Many physicists would agree with Ma's opinion that "statistical mechanics, as it is at 
present, is an ill-proportioned subject, with many successful applications but relatively 
little understanding of the basic principles" This statement can be disputed for 

classical statistical mechanics, since modern chaos theory has illuminated Boltzmann's 
"ergodic hypothesis" for closed classical systems (see, e.g., 0). However, there has been 
comparatively little work on the foundations of statistical mechanics for closed quantum 
systems. Closed quantum systems are the subject of this Letter; I demonstrate that results 
from quantum chaology (the study of quantum systems which are classically chaotic) can 
provide a sound, easily understood basis for quantum statistical mechanics. 

A closed, bounded quantum system has a discrete energy spectrum, and the state vector 
of such a system at time t can be written as 

\m) = T.Cc.e-^^-'/''\a), (1) 

a 

where \a) is an energy eigenstate with eigenvalue E^, and the coefficients Ca (normalized so 
that J2a IC'aP = 1) specify the initial state. We would like to understand how an approach 
to thermal equilibrium can be encoded in Eq. (|l|). 

To decide whether a system is in thermal equilibrium, we must measure something. Let 
us say that the quantities to be measured are represented by a set of hermitian operators Ai, 
and that our measurements can be characterized by some finite time average of the quantum 
expectation values {ilj{t)\Ai\ilj{t)). This is clearly a drastic oversimplification of any real 
experiment, but we will make the assumption that a detailed analysis of the measurement 
process is not essential to an understanding of statistical mechanics. 

Although we cannot measure them, an interesting set of objects to consider are the 
infinite time averages, 

A = }^ll^dt (^(t)IAI^(t)) 

= Y.\C^ma\A,\a). (2) 

a 

If our system approaches thermal equilibrium from a generic initial state, then we would 
expect Ai to depend on the total energy of the system, but to be independent of all other 
aspects of the initial state This is because, even though we will have some nonthermal 
behavior at the beginning (which will of course depend on the initial state), after the 
system reaches thermal equilibrium, it should stay there except for small (or large but 
rare) fiuctuations. All of these short-time effects should be washed out by the infinite time 
average. Even so, the system will have a well defined final temperature only if its total 
energy {E) = Y.a \Ca\'^Ea has a quantum uncertainty AE = {Y,a \Ca\'^E'^ — (i?)^)^/^ which 
is small: AE ^ (E). From now on we will restrict our attention to initial states which 
satisfy this condition. Returning to Eq. (|^), we can see that Ai will be independent of 
the detailed pattern of the C^s if and only if (al^jla) is a smooth function of Ea, with 
negligible variation over the relatively small range (E) ± AE. Whether or not this is a 
reasonable condition on the Aj's will be taken up later. 

For now, let us assume (Q;|y4j|a) is indeed a smooth function of E^, and write it as 
(a|y4j|a) = ^i{Ea). Thus we have 
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A = Um[^ + 0{^E/{E))]. (3) 

We can gain some understanding of the function if we consider the thermal expectation 

value 

a 

= ZiT)-' / dE n{E) e-^/^ , (4) 

Jo 

where T is the temperature (and we have set Boltzmann's constant to one), Z{T) = 
J2a e~^°'/'^ is the partition function, and n{E) = J2a ~ ^a) is the density of states. For 
a system with N ^ 1 degrees of freedom, and when integrated against a smooth function of 
E, the exact n{E) can be replaced (with negligible error) by a smooth function of the form 
n{E) ~ exp[A^s(-E/A^)], where s{E/N) can be identified as the thermodynamic entropy per 
degree of freedom. We then have 



{Ah 



J^dE e^[^(^/^)-^/^^l $,(E) 
dE e^[^(^/^)-^/^^] 

MU)[l + 0{N-'/')], (5) 



where the second line follows from the steepest-descent approximation, and U{T) = 
T^Z'{T)/Z{T) is the thermodynamic internal energy, given by the solution of the steepest- 
descent condition s'{U/N) = 1/T. Comparing Eqs. (^ and (^), and continuing to assume 
both AE <^ (E) and ^ 1, we obtain a satisfying result: the infinite time average Ai 
is equal to the thermal average (Aj)^, with the temperature T related to the quantum 
expectation value of the energy {E) by U (T) = {E) . 

Furthermore, since (al^jla) = ^i{Ea), Eq. (^ implies that (alAjlo;) = (Aj)^^, where 
Tq is given by the solution of U{Ta) = E^. Thus, if Ai is to be independent of the initial 
state, then individual energy eigenstates must have thermal properties. This phenomenon 
was previously discovered (via different arguments) in the special case of the hard-sphere 
gas, where it was dubbed eigenstate thermalization 0. 

Now let us come back to consider 

{m\A\m) = j:c:cpe^^^--''^^'/\a\Am (6) 

without any time averaging at all. We expect that {ilj{t)\Ai\ilj{t)) will start out at t = with 
whatever value is picked out by the initial state, and then evolve to its thermal expectation 
value {Ai)T after some characteristic relaxation time. Examining Eq. (^), we see that this 
can only happen if (a|y4j|/5) is generically much smaller than (a|y4j|a) when (3 ^ a. In this 
case, the (3 ^ a terms in Eq. (|^) will usually make a negligible contribution, and then we will 
find {'ilj{t)\Ai\'ip{t)) = {Ai)T, just as we found Ai = {Ai)T- However, if the magnitudes and 
phases of the C^'s are carefully chosen, then we can "line up" the {a\Ai\P)'s so as to get any 
value of {'ijj{0)\Ai\'ip{0)) that we might want. Afterward, however, as we see in Eq. (H), the 
phases will change in the usual manner; the carefully contrived coherence among the various 
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{a\Ai\P)^s will eventually be destroyed, and we will again find {■ip{t)\Ai\ip(t)) = (Aj)^. The 
outstanding problem, not to be solved here, is to understand this process more quantitatively 
for specific systems. This will entail characterizing classes of physically relevant initial 
states by some sort of statistical relationship between the C^'s and the (a|Aj|/3)'s, and then 
demonstrating that this relationship results in the appropriate phenomenological equations 
(e.g., the Boltzmann equation for a dilute gas). This program is obviously an ambitious one. 

We can, however, use our present results to resolve a paradox ^ which arises in the 
conventional approach to quantum statistical mechanics. Suppose we have a Carnot engine 
1^: a system with an adjustable parameter V (such as its volume, or an external electric or 
magnetic field) which has been brought to thermal equilibrium by placing it in contact with 
a large heat bath at temperature Ti. This means that we should describe the system with 
a "thermal" density matrix that is diagonal in the energy eigenstate basis, with diagonal 
entries given by the Boltzmann weights Z{Ti, Vi)~^e~^°'^^^'^^'^'^ , where Vi is the value of the 
parameter V. Now remove the system from contact with the bath, and change V very 
slowly from Vi to V2. According to thermodynamics, this is an "adiabatic" process, and the 
system remains in thermal equilibrium at a slowly varying temperature. At the end of the 
process, conventional statistical mechanics would assign the system a new thermal density 
matrix, diagonal in the energy eigenstate basis with diagonal entries Z(T2, V2)~^e~^"^'^'^''^'^'^ , 
where T2 is the final temperature. However, according to quantum mechanics, a sufficiently 
slow (or "adiabatic") process preserves the probability for the system to be in any given 
energy eigenstate (although the eigenstates themselves change with V^); thus the system does 
not acquire the new Boltzmann weights appropriate to its new volume and temperature, 
but instead retains the original ones. How, then, can we understand it to be in thermal 
equilibrium? 

The conventional response is to declare the question irrelevant, either because in practice 
we cannot change V slowly enough to invoke quantum mechanical adiabaticity |l| , or because 
we can only "essentially" but not "completely" isolate a system from its environment [^]. 

Our answer is different, and has already been presented: a system which exhibits 
eigenstate thermalization can come to thermal equilibrium for any initial state with a 
relatively small energy uncertainty. The discussion was for pure states, but it is trivially 
generalized to include mixed states: simply replace C*C^ by pa/3 in Eq. (|), where pa/3 is 
the initial density matrix in the energy eigenstate basis. It follows immediately from our 
previous considerations that if pap is diagonal, and if AE / (E) is small, then the system 
is automatically in thermal equilibrium. In the case at hand, the initial thermal density 
matrix (which according to quantum mechanics is unchanged during the adiabatic process) 
certainly meets these two conditions: it is explicitly diagonal, and for a system with ^ 1 
degrees of freedom we know (from conventional statistical mechanics) that AE ~ N~^^'^{E). 
However, the energy eigenstates do change with V, and so their thermal properties reflect the 
new volume V2. Meanwhile, the expectation value of the energy changes because the energy 
eigenvalues do, and this results in a new temperature. All together, this means that, after 
the slow change in V has been completed, the quantum expectation value of each Ai will be 
equal to the thermal expectation value that is predicted (for volume V2 and temperature T2) 
by conventional statistical mechanics, even though each energy eigenstate retains its original 
Boltzmann weight. In short, the paradox is resolved by eigenstate thermalization. 

All of our results depend crucially on two assumptions: that (a| Ajlo;) is a smooth function 
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of Ea, and that |(a|y4i|/3)| <^ |(a|y4j|a)| when (3 ^ a. Under what circumstances can we 
expect these assumptions to hold? A comprehensive answer is of course out of reach, since 
the only restrictions we have placed on the systems under consideration are that they are 
bounded (so that the energy spectrum is discrete) and that they have many degrees of 
freedom. Nevertheless, it is possible to identify a broad class of systems for which the 
validity of our assumptions can be demonstrated, in some cases with full mathematical 
rigor: many-body quantum systems which are classically chaotic. 

For such a system, one can make a semiclassical calculation of the matrix elements (in 
the energy eigenstate basis) of any operator Ai{p, q) which does not depend explicitly on Ti, 
and is a smooth function of the canonical coordinates and momenta. The diagonal matrix 
elements are given by |9|-|l5| 



{a\A,\a) = {A,}EA^ + 0{h)], (7) 
where we have defined the classical, microcanonical average 

_ jd'^pd^q5{E-H{p,q))A,{p,q) 
^ '^"^ JdNpd^q6{E-H{p,q)) ' 

and H{p, q) is the classical hamiltonian. Eq. (|^) has been proven rigorously when the chaotic 
system consists of geodesic motion of a particle on a compact iV-dimensional manifold with 
everywhere negative curvature P, p!4| , p!5| ; the 0{h) term vanishes as E^ oo. Clearly this 



means that (alAjla) is a smooth function of Ea for sufficiently large values of E^- We 
can, however, go further; we will see that, for any particular eigenstate |a), the expected 
deviation A(a|Aj|a) of (a|y4j|a) from the estimate of Eq. (|^) is exponentially small in N, 
and hence numerically negligible when N ^ 1. First, though, we turn our attention to the 
off-diagonal matrix elements (a|y4j|/3). 

We begin with the observation that, for a quantized chaotic system, the numerical values 
of the off-diagonal matrix elements (a|y4j|/3) should be smoothly distributed when E^ and 
Ei3 are in fixed, narrow, nonover lapping ranges. The idea here is that, since there are 
no classically conserved or approximately conserved quantities other than the total energy, 
there should be no selection rules which would make some matrix elements much larger than 



others [ll6|-[18[]. A smooth distribution for the {a\Ai\P)^s also follows from Berry's conjecture 
that the energy eigenf unctions behave like gaussian random variables p!B|JTP| -pT| . Next we 
note that if Ai is to be a useful observable, its quantum uncertainty AAi cannot be too 
large. For an energy eigenstate AAi is given by [|r^ 

{AAif = |a) - {a\Ai\af 

= EI(«I^^I/5)I' 

^niEa)Ea\{a\A,\f3)\\ (9) 

In the second line, we used completeness, and in the third, crudely approximated the sum; 
in the third line, |(a|y4i|/3)p stands for a typical value of the squared off-diagonal matrix 
element when Ep is close to E^- The key point is that n(i?) is exponentially large in N, 
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while AAi cannot be very much larger than (a|y4j|Q;). Thus, the typical value of 
must be exponentially small in A^. 

This argument can be considerably refined [p2| , p3[ ; a semiclassical estimate of the typical 
value of is 



2 Si{uj,E) 



\{a\Am' = ^^[l + 0{n)], (10) 

where u = {Ea — Eji)/h and E = ^{Ea + Ep), and we have defined the power spectrum 
Si{u!, E) of the classical observable Ai on the phase-space surface of energy E, 

1 r+oo 

S,{u;,E) = — / dte'-' {{A{t)AmE - {Ami) ■ (H) 

Again, the key point is that Si{uj,E) can refiect only the intrinsic dependence of Ai 
(e.g., extensive or intensive), while n{E) grows exponentially with A^; thus the off-diagonal 
matrix elements must be exponentially small in A^. This conclusion must hold even after 
all quantum corrections to Eq. ( p!OD are included. (Also note that the factor of 1/h on the 
right-hand side of Eq. (^) does not mean that the off-diagonal matrix elements become 
large when h is small, since l/n{E) is proportional to .) 

Now we return to our discussion of the diagonal matrix elements (a|y4j|a). For any 
particular eigenstate \a), the expected deviation A(a|y4i|a) of (alAjla) from the estimate 
of Eq. (0) is given by the u; ^ limit of Eq. (^) ||; that is. 

The previous discussion of A^ dependence applies, and so we see that A(a|y4j|a) is 
exponentially small in A^. Once again, this conclusion must hold even when all quantum 
corrections to Eqs. (^ and {\L2^ are included. Thus, for a quantized chaotic system with 
many degrees of freedom, we expect (alAjla) to be a smooth function of Ea for any operator 
Ai{p,q) which is a smooth function of p and q, and we also expect |(a|74j|/3)| -C |(a|74j|a)| 
when (3 ^ a. These are precisely the conditions we need in order to understand the 
emergence of statistical mechanics from the underlying quantum dynamics. 

To conclude, we have seen that a closed quantum system with many degrees of freedom 
can approach thermal equilibrium from a generic (but definite) initial state if the observables 
Ai which are measured to check for thermal equilibrium satisfy two conditions. First, the 
expectation value (al^jla), where \a) is an energy eigenstate, must be a smooth function of 
the energy eigenvalue E^- Second, the off-diagonal matrix elements (a|y4j|/3) must be much 
smaller than the corresponding diagonal ones. If the first condition is met, then (Q;|Aj|a) 
turns out to be equal to the thermal expectation value of Ai at a temperature Tq, given by 
U{Ta) = Ea, where U{T) is the thermodynamic internal energy of the system (computed in 
the usual way) as a function of temperature; this is eigenstate thermalization. If the second 
condition is also met, then we have a basis for understanding how a system which is prepared 
in a nonequilibrium state can evolve to thermal equilibrium via hamiltonian time evolution. 
Both conditions are expected to be satisfied by all simple observables whenever the quantum 
system is obtained by quantization of a many-body classical system which is chaotic. Of 
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course, the two conditions could also be satisfied even if the quantum hamiltonian does not 
have a classical limit (e.g., a lattice spin system, perhaps with the requirement that the Aj's 
are local or bilocal). An interesting open question is whether or not quantum systems of 
this type have any other properties in common (such as eigenvalue spacing statistics p4| ) 
which would serve to identify them more easily. 
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